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ABSTRACT 



The illumination condition of an artificial earth satellite 
is investigated in relation to the geographic location of the 
sub-satellite point. A stochastic problem is formulated by 
postulating a *’random*^ satellite, and the solution to this 
problem is indicated analytically. A statistical solution is 
obtained from a digital computer simulation program, together 
with a sensitivity analysis of parametric variations. 
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PREFACE 



The general purpose of this thesis is to investigate certain aspects 
of the illumination of artificial earth satellites by the sun j in relation 
to the track of the sub- satellite point as measured in a geographic 
coordinate system. Imposed upon this study is a stochastic model 
of the earth- satellite system. 

An analytic approach to the basic problem was made , followed 
by a computer simulation program for a statistical solution. The 
basic problem is set forth in the Introduction. 

I wish to express my appreciation to the following persons for 
their contributions to this study: 

Mr. E. J. Leinfelder , Chance- Vought Corp. , Astronautics 
Division, who first posed the problem to me, and who provided the 
stimulus and inspiration leading to this study; 

Professor W. P. Cunningham, U. S. Naval Postgraduate School, 
who indicated the method of formulating the problem in a manner 
which led to the desired results; 

Professor R. M. Thatcher and Miss Pat Hoang, U. S. Naval 
Postgraduate School, who assisted me in no small measure in the design, 
writing and de -bugging of the computer simulation program; and 
Professor R. R. Read, U. S. Naval Postgraduate School, who 
detected initial errors in the chapter on probability formulation. 
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Symbol 

a 



dl 

h 



W 

I’t 

I’o 



r 

r 

s 
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A 

C 

D 



Definition 

semi-major axis of orbital ellipse 

apparent change in longitude of sub- satellite 
point from initial reference time to time t 

altitude of satellite above surface of earth 

inclination of satellite orbital plane to 
equatorial plane 

inclination of equatorial plane to ecliptic plane 

longitude of designated geographical points 

longitude of sub- satellite point at time t 

longitude of ascending node 

unit vector from center of earth to satellite 

projection of r onto equatorial plane 

specified angular distance of sub- satellite point 
from designated geographical points 

instantaneous time after initial reference time 

angular rotational velocity of earth about polar 
axis in radians per unit time 

angular velocity of satellite in radians per unit 
time 

angle between E and r vectors 

angle at satellite subtended by the center of the 
earth and the surface of the earth 

angle between vernal equinox and ascending node 
of satellite 
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Definition 



Symbol 



E 

F 

G 

G' 

H 

L 

L' 

M 

M' 

N 

Re 

T 



unit vector fronn center of earth to sun 

angle between ascending node and projection of r 
onto equatorial plane 

angle between ascending node and r vector 
universal gravitational constant 
angle between vernal equinox and E vector 
latitude of designated geographical points 
angle between r and r' 

event that desired passage occurs during first orbit 
mass of the earth 

event tha,t desired illumination occurs at time of 
closest approach to designated geographical point 

radius of the earth 

orbital period of satellite 
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INTRODUCTION 



It is proposed to investigate some aspects concerning the 
illumination of an artificial earth satellite by the sun in relation to 
the track of the sub- satellite point as measured in a geographic 
coordinate system. But before stating the complete problem, some 
preliminary discussion is necessary. 

Consider a satellite in a circular orbit around the earth. Its 
position may be completely specified at any time (ignoring 
perturbations) if its altitude, orbital inclination to the equatorial 
plane, and the longitude of its ascending node at a particular time is 
known. The longitude of the ascending node may be defined as the 
longitude at which the sub- satellite point crosses the equator entering 
the northern hemisphere. Under these circumstances the track of 
the sub- satellite point on the earth *s surface may be computed for any 
set of sequential orbits. 

Due to the rotation of the earth about its polar axis, the longitude 
of the ascending node is displaced westward on each orbit by an 
amount equal to the rotational velocity of the earth times the orbital 
period of the satellite. This displacement will be denoted w^T , where 
Wg is the angular rotational velocity of the earth about its polar axis, 
and T is the orbital period of the satellite. 

If it is desired to determine if the track of the sub- satellite 
point will pass within a given distance of a given point on the earth *s 
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surface, then under the stated conditions the problem is completely 
deterministic. A simple but somewhat tedious calculation will 
ascertain the number of orbits required, the exact time, and the 
actual distance from the given point at the time 6f closest approach. 

Now suppose that, instead of being completely specified, these 
parameters of altitude, inclination, and longitude of ascending node 
were random variables. Conceivably this situation could arise where 
a satellite has been ^'lost^’ to surface detection units for an extended 
time period. Similarly, this proposition might be required in a 
feasibility study of a ground visual detection station. In any event, 
under these circumstances the problem of ascertaining if the satel- 
lite track passes within a given distance of a given point is no longer 
deterministic but stochastic. The solution must be given as a 
probability of such passage after a given number cf orbits, and will 
be a function of these three random variables. 

To simplify the problem considerably, let the variables of 
altitude and inclination take on a range of discrete values. Then seek 
a solution for the random variable of longitude of ascending node for 
each pair of values for altitude and inclination. The longitudes of 
the second and all subsequent ascending nodes are simple functions 
of the first nodal longitude, since specifying the altitude specifies 
orbital period, hence specifies the longitude displacement, w^,T. 
Then the probability of passage is a function of the random variable 



2 










■"M 






1 1 : 4 



# 




• II I 



1 






of longitude of first ascending node, and of the orbital period, T. 

The second part of the problem is to determine if the satellite 
is also illuminated by the sun. It appears obvious that the satellite 
need not be over the sunlit side of the earth to be illuminated. In 
fact, the only requirement for illumination is that the satellite be 
outside the shadow cone of the earth. For satellites at sufficient 
altitudes this may not occur until the satellite is well into the ’’back 
side’* of its orbit. 

With these considerations in mind, the main problem may now 
be stated: 

Given a satellite at an altitude h, with an orbital inclination i, 
and two geographical points on the earth’s surface on the same 
meridian of longitude W, and equidistant from the equator at latitudes 
L North and South: 

Find the joint cumulative probability that the satellite will pass 
within a given angular distance s of either point and that the satellite 
will simultaneously be illuminated by the sun, the result to be stated 
as a cumulative probability distribution function in terms of successive 
orbital periods after initial reference time (arbitrarily set equal to 
zero). 

Some previous studies of satellite illumination have been made, 

made, but no studies have been found which treat the problem from a 

1,2 

stochastic viewpoint. 



3 









< 



I — l>W1— 



I 



<ir; 



«i«pf ^ 

-»HU‘ 

4 * 






This thesis first develops the necessary mathematical models 
for the orbital mechanics, then sets forth the precise probability 
formulation statements. The statistical solutions are then derived 
from a digital computer simulation program which is discussed in 
the appendices. 

While the design and writing of the computer simulation program 
proved of great interest to the author, and provided an opportunity 
to apply recently acquired knowledge of Monte Carlo techniques, it 
must be emphasized that the computer program is not an entity in 
itself, but merely a useful tool in arriving at a solution to the stated 
problem. 



Leinfelder , E.J. , and Thrasher, J.R. , Orbit-Earth Shadow 
Relations, AST/EOR- B277 , Chance-Vought Corp. , Dallas, Texas, 
1960. CONFIDENTIAL 

^Koelle , H.H. , Editor, Handbook of A stronautical Engineering , 
McGraw-Hill Book Co. , New York, 1961 . 
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CHAPTER I 



ORBITAL MECHANICS 



Initial Assumptions 

!• The earth *s solar orbit and the satellite orbit have zero 
eccentricity. 

2* The earth is a perfect sphere with total mass concentrated 
at the center. 

3. The earth *s equatorial plane is inclined to the ecliptic 
plane at a constant angle i' (= 23^ 26* 59*0 • 

4. The earth*s solar umbra is a circular cylinder (valid for 
values of h less than three times the earth's radius). 

5. At initial reference time the longitude of the ascending node 
of the satellite may be considered a random variable , uniformly 
distributed on the interval ^0,27rj . At this time the satellite will 

be assumed to be located at the ascending node. 

6. All other perturbations of a satellite orbit will be ignored. 
This includes atmospheric drag and solar radiation pressure. 

Discussion of Assumptions 

The effect of these assumptions is to confine the problem to 
one involving only central forces. This also implies that precession 
of the satellite orbital plane will not be included. In order that these 
omissions do not radically affect the solution, the probability dis- 
tribution function sought will only be examined for a period of 16 orbits. 
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This will be on the order of 24 hours for the range of altitudes con- 
sidered. As will be seen subsequently, two more random variables 
must be introduced; one to account for the space orientation of the 
satellite orbital plane, and another to account for seasonal variations. 

Method of Attack 

A mathematical model will be established representing the motion 
of the satellite and the sun with reference to a geocentric coordinate 
system, followed by a probability formulation. 

Mathematical Model 

Let E represent the vector from the center of the earth to the 
sun, and let r represent the vector from the center of the earth to 
the satellite. The initial assumptions will ensure that these vectors 
are of constant magnitude and possess constant rotational velocities 
about the center of the earth. 

Using Figure 1“ 1 , find the direction cosines of E in the rectan- 
gular geocentric coordinate system depicted, where the X-axis is 
along the vernal equinox, the Z axis is directed through the north 
celestial pole, and the Y-axis perpendicular to the other two axes in 
a conventional manner. The origin is, of course, the center of the 
earth. 

During a yearly cycle the angle ^ takes on all values in the 
interval I 0,271 . The values of x, y, and z are easily determined as: 
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Fig. I'i; Representation of 
Earth -Sun vector E 
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X = E cos H 



\ 



y = E cos i * sin H 



(M) 



z “ E sin i * sin H 



Hence the direction cosines are : (cos H, cos i* sin H, sin i* sinH). 

From Figure 1-2 find the direction cosines of r in the same 
coordinate system. The angle D may be thought of as the right 
ascension of the ascending node. The angle G_ represents the 
angular distance of r from the ascending node. Obviously the angles 
D and G take on all value in j^O >27tj . 

From Figure 1-2 the angle F is the projection of G onto the 
equatorial plane. Let r* be the projection of r onto the same X-Y 
plane. Then r* - r cos (angle between r and r*). 

Using Figure 1-3 the right spherical triangle with sides G, F and 
L * reveals that L ’ is the angle between r and r*. Then 



sin L* “ sin i sin G 



(L2) 



and 




Also from the right spherical triangle in Figure 1-3: 



tan F = cos i tan G 



(M) 



hence 



F - tan ^ (cos i tan G) 



(1-5) 
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Fig. Representation of earth 
Satellite vector F 
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Ftg. I“3 : Simpl[fi.ed view of 

Earth " Satellite vector r 
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Referring again to Figure 1-2, 



X “ r * cos (D+F) 

- r cos L ’ cos (D+F) 

= r (1 - sin^ i sin^ G) ^ cos 

y = r * sin (D+F) 

/I . 2 . . 2 

“ r (1 - sin 1 sin G) sin 

z = r sin L * 

- r sin i sin G 



-1 

D + tan (cos i tan G) 

( 1 ^) 



D + tan^^(cos i tan G) 
(1-7) 



(1-8) 



Using Equations (1-6), (1-7) and (1-8), the direction cosines of r are: 



. 2 . . 2 

(1 - sin 1 sin G) cos 

2 9 

(1 - sin i sin^ G) sin 



D + tan" ^ (cos i ta,n G) 

-1 



D + tan (cos i tan G) 
sin i sin G ^ • (1-9) 

Referring to Figure 1-4 it is noted that the angle G describes 
the angular position of the satellite (or the sub- satellite point) in the 
orbital plane from the ascending node at time t. Hence 

G = Wgt (i-10) 

where 



Wg = angular velocity of satellite in radians /unit time 
t “ instantaneous time after reference time. 

As the satellite proceeds in its orbit to time t, the earth also rotates 
about its axis with an angular velocity . Then dl is the apparent 
change in the longitude of the sub- satellite point in time t. Clearly, 
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Fig. 1“^: Location of sub~ satellite 
point at time t 
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F = Wgt + dl 



Then the longitude of the sub- satellite point at time t is; 



IV = 1' + dl 

t o 



I'o + F - Wgt 



1*^ + tan ^ (cos i tan G) — w^t 



( 1 - 12 ) 



It is more convenient at this point to develop an analytic model 
which will indicate whether or not the satellite is illuminated by the 
sun. Figure 1-5 depicts the plane containing the E and r vectors. 
For simplicity, let both E and r be unit vectors. 

As the satellite orbits about the earth, the angle ^ (between 
the E and r vectors) can take on all values in the interval j^0,27xj 
However, since the representation in Figure 1-5 is symmetric about 
the E vector, the problem requires an interest in the magnitude of 
A only in the interval 



If cos A is non-negative, the satellite v/ill be illuminated by the sun. 

If cos A is negative, the satellite may or may not be illuminated. 

Let B =A-7X/2as shown in Figure 1-5. Let ^ be the angle at 
the satellite subtended by the center of the earth and the surface of the 
earth. Then 




Clearly , 



E • r 



cos A 



(1-13) 




(1-14) 
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Fig. I"5 : Plane containing E and r vectors 
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a constant for circular orbits, where is the radius of the earth. 
Then if B + C is greater than TT/Z, the satellite will be in the 
shadow cylinder of the earth. Hence if the satellite is illuminated, 
then 



or 



cos (E • r) 



7t 

2 



+ sin 



R 



Re + h 




(1-15) 



cos 




r) 



. -1 
+ Sin 




= IT 



( 1 - 16 ) 



If it is not illunninated , then 



cos ^(E • r) + sin" ^ 



R 



+ h 
e 



> TT 



(1-17) 



Using the results of Equations (1-1) and 1-9), then: 
cos A = E • r = (1 - sin^ i sin^ G) ^ X 



{cos H cos [ D + tan ^(cos i tan G) j + 
cos i ' sin H sin D + tan” ^ (cos i tan G)j | 



+ sin i * sin i sin H sin G 



(1-18) 



The next stage of the problem is to determine if the satellite 
track passes within a given angular distance s of either of two 
geographical points with coordinates (L, W) and W). With this 

concept the satellite track may be considered to sweep out a ’’swath*’ 
over the earth’s surface, of width 2s, Denote 1*^ as the longitude 
of the ascending node. 
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For illustrative purposes it is advantageous to show a schematic 
plot of the satellite track across the surface of the earth, with 
special emphasis upon its relationship to the designated geographical 
points. For clarity, a Mercator projection will be used. 



For cases where 



- s > 



then the situation in 



Figure 1-6 will exist. Clearly, if the value of W lies in one of the 
intervals [ ^ j » ^ + l] J “ 3, 5, 7, then the desired passage 



will occur. Indeed this is a necessary condition 

< 



For the case where 



< 

S 



+ s , then Figure 



1-7 will apply. Here, the desired passage will occur if W lies in 
one of the intervals ^Ij , ^-flj > j “ 1, 3. This too is a 

necessary condition. 

Using these conditions to determine if passage occurs, the time 
t of closest approach to the geographical point may be evaluated. 
Using this value of t to compute cos A in Equation (1-18), then the 
conditions contained in Equation (1-16) are applied to determine if 
the desired illumination occurs. 

It is now necessary to develop expressions for the orbital period 
and angular velocities for the satellite and the earth. It can be shown 



that the period of a satellite is: 
T = 2 7T 



(1-19) 
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t 







where 



M 

In mks units , 



a = semi-major axis of the orbital ellipse 



= + h, for circular orbits 



= universal gravitational constant 
= mass of the earth. 



Re 


= 6.371 


X 


10^ 


/ 

G 


= 6 . 668 


X 


10"^^ 


/ 

M 


= 5.98 


X 


10^^ ! 



Then 



= 2 71 (6.371 X 10^ meters + h) 



Ml 



6.668 X 10"^^ nt.-m^ X 5.98X10 



24 



kg. 



kg' 



3/2 



= 2 -n (6.371 X 10^ m. + h) 

(3.99 X 10^^ /sec^) 

= 3. 148 X lO"”^ (6.371 X 10^ + h) sec 



= 5.246 X (6.371 + h) 
>6 



3/2 



1/2 



min. 



( 1 - 20 ) 



where h is given in units of 10 meters. 
Since w _T = 2 tt , 

o 

Wg = 2 7T /T 

and 



Wg = 2 TC /Tg 



( 1 - 21 ) 



where is the rotational period of the earth and 
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( 2 TC / 1440 ) min . 

4.363 X 10~^ min. (1-22) 
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CHAPTER II 



PROBABIEITY FORMULATION 

The next stage leading to a solution of this problem is the 

formulation of a probability statement. The probability sought is 

the probability that illumination occurs simultaneously with the 

tVi 

desired passage by the orbit, for n = 1, 2, . . . , 16. 

Looking at the first orbit, the longitude system may, without any 
loss of generality, be arbitrarily defined such that the longitude of the 
ascending node, 1*^ , is equal to zero, and all other longitudes are 
measured eastward from zero through a total angle of 2 TC . 

Let M be the event that the desired passage occurs sometime 
during the first orbit, and let N be the event that the desired 
illumination occurs at the particular time t when the satellite makes 
its closest approach to the geographical point. Note that M is a 
time-independent event, while N is strictly dependent upon the value 
of t. The probability sought is P j^MNj . 

Studying each event separately, it is clear that 

P [ m] = P [ Ij = W = ]j + i] (2-1) 

where j = 1 and 3, or else j = 1, 3, 5,7, depending upon the values 
of i , s and L. 

This equation can be written in more meaningful language if 
P Mjj is defined as the probability that W lies in ^ Ij , ^ + l] 
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Clearly, the events Mj^, M 3 , . . . , are mutually exclusive, 
hence 

P [m] = ^ P [Mj] (2-2) 

odd j 

Also, from Equation (1-16), 

p[n]=p[a + c =n-] =p[a in . c] (2-3) 

Hence the random variable associated with event M is W, and 
it appears intuitively correct that A is the random variable associated 
with event N. Then the joint probability sought is 

P [ ] (2-4) 

odd j 

To verify that A is a random variable, note that Equation (1-18) 
shows that A is a function of i, G, D and H. But for a particular 
time t, the values of i and G are constant. Then 

is simply the probability that D and H take on a range of values such 
that A = TT — C , which implies that D and H are random variables 
with associated probability density functions, f£)(d) and f jj(h) , both 
defined on the interval 0 , 2 tT ] . Care should be taken to avoid 
confusion with h, the satellite altitude; however, in most situations 
the meaning will be clear from the context. 

Using the basic assumptions and the physical interpretations of 
D and H, it appears intuitively correct to state that they are inde- 
pendent, uniformly distributed random variables such that 



[ MN ] = y 
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HI 



£„(d) 



fH(h) 



for 0 = d ^ Znl 



elsewhere 



for 0 = h s 2n\ 



elsewhere 



(2“5) 



( 2 - 6 ) 



It quickly follows that A = A{G,H) is a random variable with 
an associated probability density function, fy^(a), defined over the 
interval ^ • 

While it appears unlikely that the mathematical expression for 
f^(a) could be derived quickly by analytic methods, the formulation 
of the problem will continue as if were completely described. 

Before proceeding further with Equation (2-4) it is necessary to 
develop the form of the probability density function for W. Since the 
only quantity of interest is W - Vq , and since one of the basic initial 
assumptions was that 1*^, was uniformly distributed, it follows that 
W is also uniformly distributed. Arbitrarily setting equal to 
zero, as before, merely ensures that W is defined over the interval 
[ 0,2 71 ]. 

Then 

0 elsewhere 



l/2 7t for 0 ^ w ^ 2 tt 



(2-7) 



Now by the rules for conditional probabilities 



[hlN] 



= P 



N 



Mj ] p[mJ 



(2-8) 
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and 



[mn] 



odd j 



[ N Mj ] . p[Mj] 



(2-9) 



Assume now that, associated with P j^N Mj j , there e:d.sts a 
conditional probability density function 1^(9- | '^) » defined over the 
interval , Tt j . Then Equation (2-9) may be written as 

r 1 Vv I 

P j^MN J ~ y I I 



odd j !• 



1 

2 n 



lj + 1 TC-C 

^ ( 2 - 10 ) 
odd j ]j 0 



using the results of Equation (2“7). 

The final probability to be calculated is the joint cumulative 

probability that the desired events M and N will occur simultaneously 

th 

at least once by the " orbit, for n - 1, 2, . . . , 16. To accomplish 
this requires some new notation. 

Let represent the desired probability of at least one such 
joint occurrence by the orbit. Let P(l) be the probability that 
the desired event occurs on the first orbit; clearly P(i) is the solution 
to Equation (2-10). Let Pfu^jn+l) be the joint probability of no such 

3 1 

occurrence on the first n orbits and one occurrence on the (nf 1) 

orbit. 
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Then 



Pj = P(l) 

P 2 =Pi+P(l',2) 

P 3 =P 2 + P(2',3) 
and, in general, 

^n+1 “ P(n',n+1) , (2-11) 

This equation describes symbolically the procedure for 
calculating for n = 1, 2 , . • • , 16. The problem of calculating 

the probabilities for orbits subsequent to the first can be considerably 
alleviated by utilizing a characteristic of the satellite tracks. 

From Equation (1-12) it is easily shown that the longitude of the 
ascending node for the second orbit is 

i.e. , r»p has been displaced westward of 1 *^ by the amount w^T. 

But this is equivalent to retaining 1*^ as the longitude of the ascending 
node for subsequent orbits and displacing W eastward a distance w^T. 
In other words, Equation (2-10) may be used for the n^^ orbit, pro- 
viding the longitude of the geographical points is set equal to 

(W)^ “ W + (n - 1) w^T for n = 1,2, . . . , 16 (2-13) 

Nonetheless , there remains unsolved the problems of determining 
the limits of integration (with respect to w) and discovering the 
analytic expression for Equation (2-10). These problems 
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will be avoided by abandoning the analytic approach in favor 
of a computer simulation program using Monte Carlo method 
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CHAPTER III 



SIMULATION PROBLEM 

It is expedient first to examine the form of the desired solution 
in order to assess the proper method of arriving at it. 

The cumulative probability distribution function sought is not a 
continuous function, but a monotonic step function, with certain 
probabilities associated with each point corresponding to the end of 
an orbital period. In fact, the desired result will be a distribution 
function truncated at 16 orbital periods. 

Now assume for the moment that the problem can be set up for 
simulation, where a suitable number, k, of samples can be obtained 
from a pseudo-random number generator. Let 

Xj = number of times that passage occurs on the 1®^ orbit 
xj = number of times that passage occurs on the i^^ orbit 

and not on the i— 1 previous orbits , for i == 2 ,3 , ... 16. 

s t 

yj = number of times that illumination occurs on the 1 orbit, 
given that simultaneous passage occurs on the first orbit. 
y£ = number of times that illumination occurs on the i orbit, 
given simultaneous passage on the same orbit, and that 
such joint event did not occur on the i—1 previous orbits , 
for i == 2 ,3 ,. . . 16. 

Here , ''passage*' implies that the satellite track passes within the 
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specified angular distance s of either of the two geographical points. 
With these precise definitions it is possible to assert that: 
e[Pi] = E [ P(l)] = 

“TT 

E [Pj] ? yj + E [P(l'.2)] ? yj + y^ 

"X* ____ 

and, in general, 

IJ, ''i 

If P^(p) is defined to be the cumulative probability of passage 

(independent of illumination) , then, in similar manner, 

E [Ppjp)] ^ iV X. (3„2) 

kZ_.i=i 

Now define u ,-2 and v-.- such that 
ij ij 



u. 



ij 



V. 



ij 




0 

= / 1 



if the desired passage occurs for n 

4-lU 

the first time on the i “ orbit, 
th 

for the j sample random number > 

value 

Otherwise 

If the desired illumination and n 



(3-3) 




passage occurs for the first time 
on the i^^ orbit , for the sam- 
ple random number value. 
Otherwise 



(3-4) 
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where j = 1,2,. . . ,k. Note that v-j 
Then k 




1 does not require that ujj = 1. 

(3-5) 



and 









(3-6) 



Hence n 

i=i 

and 

E [Pn(P)] = i 

Equations (3-7) and (3-8) can be considered as the basis for the 
siniulation model. 

The design of the computer simulation program is described in 
some detail in Appendix A , with a description of some of the pro- 
cedures used to reduce the required computation time. The language 
used in writing this program was the CDC-1604 version of NELIAC, 
which is peculiarly well-adapted to this type of program. 

Appendix B contains a summary of the output of this program, 
along with a sensitivity analysis for the effect of single parametric 

variations of satellite altitude, orbital inclination, latitude of geo- 
graphical points, and angular distance s. 
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APPENDIX A 



DESIGN OF A COMPUTER SIMULATION PROGRAM 



The primary concept of the program is an artificial-time- 
controlled sequence of computations, based upon Equation (L16) and 
a modification of Figures 1-6 and 1-7. The satellite is figuratively 
put in orbit at its ascending node, corresponding to time t = 0, and 
where the longitude 1*^ is also identically zero. A auitable time 
increment is used to activate an artificial time *'clock’' which 
determines the instantaneous latitude and longitude of the sub- 
satellite point. 

Certain short-cuts are immediately available to reduce the 
otherwise inordinate number of calculations required. 

Figure A-lis similar to Figure 1-6, except that the end points of 
the longitude intervals correspond to the instantaneous longitudes of 
the satellite track when the instantaneous latitudes equal either L - s 
or L + s. Clearly the desired passage cannot occur unless the 
instantaneous track latitude L*^ lies in the interval ^L-s, » 

which implies that if the instantaneous longitude 1*^ lies in the interval 
for j = 1, 3, 5, 7, then passage will occur if W also lies 
>vithin this same interval. 

In similar manner, the interval ih- ■i+ 3 ] for j -1,5, is defined 
for situations where L-s=i = L+s. See Figure A-2. Again, 



A-1 




A -2 




A -3 



clearly, passage will occur if 1'^ and W lie within this interval. The 
numbering of the subscript j is designed to facilitate the computer 
programming . 

From Equations (1-2) and (1-10): 
sin L' = sin i sin Wgt 



hence for L* = L-s, 



t = J. 
w 



sin"“^ 



sin (L — s) 



“ ti 



sin 1 



(A-1) 



and for L*t = L+s , 



t =J. 
w. 



sin"^ 



sin (L + s) 






sin 1 



(A -2) 



where tj corresponds to the time at which the satellite track attains 
a value of longitude equal to Ij . The remaining values of tj for j 
2 can be determined quickly by 



ts ^ 

u = 

t5 = 
^6 = 

^8 




^2 

^1 

^2 

^3 

^4 



which can be easily verified by the symmetry of the track. T is of 
course the satellite period. 

Using Equation (1-12) it follows that 
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(A -3) 




and accordingly 




(A -4) 



where Ij and are defined as in Figures A-land A-2. In these 
figures if W lies within any of the intervals depicted then passage 



the same manner as Equations (A -3) and (A^). 

Then the first test for the desired passage is determining if the 
sample value for the random variable W lies within one of the afore- 
mentioned intervals. If W satisfies this test, then the test for illumi- 
nation may be made. However, it is possible for W to be outside 
such intervals and still provide for the desired passage. If W lies 
outside of, but sufficiently close to, either end point of the depicted 
intervals, then passage may occur as in Figure A-3. 

At this point it is appropriate to introduce the approximation 
' of the true instantaneous angular distance between the sub- satellite 
point and the designated geographical point. For values of s' 
sufficiently small, the curved surface of the earth can be approximated 
by a plane. Then, from Figure A-3, 





1/2 



(A -5) 
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Fig.A'3* Mercator projection of tracl<j 

special case 



A-6 



where the cosine term is used to compensate for the ^‘shrinkage*’ 
in distance between adjacent meridians of longitude with increasing 

latitude . 



cos L* 



It should be obvious that passage cannot occur if 
> 



(W -1'^) X 



= s , hence this quantity may be used as an upper bound 
in an iterative procedure for computing successive values of s' for 
successive values of time. Due to the monotonic nature of the 
satellite track, a minimum value of s' will exist for some value of t, 
say t'<^. This value of s' is compared to s; if s' - s, then the desired 
passage will occur. Otherwise, passage cannot occur. 

If these tests are applied to the i^^ orbit and j^^ sample value of 

the random variable W, the value of u^j in Equation (3-3) is determined 

If Ujj = 1, then the illumination condition may be determined by 

application of Equation (1-16). Applying this test to the i^^ orbit and 
th 

j ” sample values of the random variables W, D and H determines 
the value of Vj^j in Equation (3-4). 

However, it may happen that the desired passage will occur with 
no illumination on an earlier orbit, but with the desired passage and 
illumination on a later orbit. In order to compute Pp(p) , it is 
necessary to record the earlier event where u^j = 1, but continuing 

the process in order to compute P^. 

Since this simulation program was only incidental to the solution 
of the main problem, the computer flow diagrams and NELIAC 
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language flow charts have not been included. However, the IBM 
binary card decks and program print -out are available from 
Professor Richard M. Thatcher, of the Operations Analysis 
Department. 
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APPENDIX B 



SUMMARY OF PROGRAM OUTPUT 

The parametric inputs to the simulation program were satellite 
altitude h, orbital inclination i, latitude of designated geographical 
points L, and angular distance from these points s. The values 
of h and s are given in units of nautical miles , but were converted 
internally to units of 10^ meters and radians respectively; the 
values of i and L are given in degrees, but were likewise converted 
internally to radians. The values of these parameters and the in- 
cremental values used in the program are given below: 
h: 100(300)1900 nautical miles 

i: 20®(8°)68° 

L: 25°(10°)65° 

s: 100(100)500 nautical miles 

The value of k (sample size) was hopefully selected to be 200 , 
in order that the approximations contained in Chapter III might be 
statistically valid, and yet require a reasonably short amount of 
computer time. Even so, about three and one-half hours actual 
operating time was consumed on the computer for a single run 
through all values of the above parameters. 

One immediate result from the output of this program was the 
discovery of a "semi-optimal" distance s = 300 nautical miles. 
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In many instances, for various values of h , i and L, the cumulative 



probabilities revealed little or no increase with values of s greater 
than 300. In a significant number of cases the rate of increase in 
probability declined sharply after s = 300; it is to be noted, however, 
that this was not always the case. Frequently the rate of increase 
in probability with increasing s was approximately linear. 
Nevertheless , it was decided to limit the tabulated values of the 
probabilities to those computed for s = 300, To have done otherwise 
would have increased the size of this appendix by a factor of five. 

Another unexpected discovery concerned certain substantial 
decreases in cumulative probabilities for h = 1900 nautical miles. 

The reduction did not occur consistently for all values of the other 
parameters, however. It is suspected that this reduction is most 
likely a statistical anomaly; unfortunately insufficient time was 
available for further repeated runs using this value of h. This so- 
called anomaly is suggested as a topic for further study, as is an 
analysis of variance for repeated runs. The computer program is 
easily adapted to the study of the effects of variation of any one 
parameter. 

Figures B-1 through B-10 give the joint probability of passage 
and illumination vs. number of orbital periods, while Figures 
B- 11 through B-15 give the independent cumulative probability of 
passage vs. number of orbital periods. 
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Figures B-1 through B-4 and B-11 through B-15 represent 
probabilities averaged for all the altitudes examined and, as such, 
show the sensitivity due to variations in inclination. The rationale 
for this averaging process lies in the fact that the variations in 
probabilities due to increasing altitude were sufficiently small to 
warrant doing so. 

Figures B-5 through B-ID , however, are given for specific 
latitudes and inclinations, and show the sensitivities due to 
variations in altitude. In these cases it will be noted that the values 
of L and i are approximately the same. In these cases the probabi- 
lities are consistently higher than in cases where i is considerably 
greater than L. Of course , when i is less than L by 3^ or more , 
no passage occurs. 

These graphs have been smoothed to some extent over the raw 
statistical output, but the general trend has been retained in all 
cases. No attempt should be made to interpolate between discrete 

I 

orbital periods. 
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